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10.  ABSTRACT  (Centlwte  on  teeetee  elde  II  neeeee.ur  "ltd  I  dent  tty  kr  Heed  mMiSor? 

A  method  of  suppressing  oce^r  surface  wave-generated  magnetic  noise  froi.  near  surface  mounted  total  field 
magnetometer  signals  is  presented.  This  technique  uses  reference  data  furnished  by  a  measurement  of  the  ocean  surface 
displacement  at  a  single  point  in  the  vicinity  ol  the  magnetometer  to  cancel  adaptively  the  magnetic  noise  generated  by 
the  surface  wave  field.  A  transfer  function  is  derived  which  relates  the  ocean  surface  displacement  measurement  to  the 
associated  magnetic  noise  field.  Experimental  data  obtained  front  a  total  field  magnetometer  mounted  7  meters  above 
the  ocean  surface  in  approximately  18  meters  of  water  display  good  correlation  with  an  ocean  surface  displacement 
measurement  made  with  a  sub-aurface-mounted  pressure  gauge.  The  transfer  function  which  relates  there  measure¬ 
ments  shows  good  agreement  with  the  theoretically  derived  transfer  function.  When  the  surface  displacement  time  . 
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PROBLEM 

Develop  a  method  of  suppressing  ocean  surface  wave-generated  magnetic  noise  from 
near  surface-mounted  total  field  magnetometer  signals. 


RESULTS 

Reference  data  furnished  by  a  measurement  of  the  ocean  surface  displacement  at  a 
single  point  in  the  vicinity  of  the  magnetometer  were  used  to  cancel  adaptively  the  magnetic 
noise  generated  by  the  surface  wave  field.  A  transfer  function  is  derived  which  relates  the 
ocean  surface  displacement  measurements  to  the  associated  magnetic  noise  field.  Environ¬ 
mental  data  obtained  from  a  total  field  magnetometer  mounted  7  meters  above  the  ocean 
surface  in  approximately  1 8  meters  of  water  displayed  good  correlation  with  an  ocean  surface 
displacement  measurement  made  with  a  subsurface-mounted  pressure  gauge.  The  transfer 
function  which  relates  these  measurements  showed  good  agreement  with  the  theoretically 
derived  transfer  function.  When  the  surface  displacement  time  series  was  used  as  the 
reference  input  to  an  adaptive  noise  canceller,  the  surface  wave-generated  noise  was 
suppressed  by  approximately  20  dB  at  the  main  swell  frequency  of  0.08  Hz.  Within  the 
frequency  band  between  0.05  Hz  and  0.3  Hz,  the  degree  of  cancellation  varied  between 
10  d3  and  20  dB. 


INTRODUCTION 

Onr  approach  to  the  problem  of  detection  and  localization  of  submerged  objects  is  to 
take  advantage  of  the  magnetic  fields  generated  by  iron-bearing  material  contained  within  the 
object  of  interest.  This  magnetic  signature  may  be  obtained  with  one  of  several  types  of 
magnetometers  whose  sensitivities  range  from  about  1  nT*  to  better  than  10“*  nT. 

To  take  full  advantage  of  the  inherent  sensitivity  of  one  of  the  better  instruments, 
one  must  be  able  to  cope  with  a  number  of  noise  sources,  any  one  of  which  may  mask  the 
additional  capability  of  a  more  sophisticated  magnetometer. 

For  a  stationary  sensor,  two  noise  sources  predominate.  The  first  of  these 
components,  commonly  called  ‘geomagnetic  noise',  has  its  origins  in  a  complicated 
interaction  of  the  solar  plasma  with  the  earth's  magnetic  field.  Charged  particles  become 
trapped  and  circulate  in  a  variety  of  trajectories,  causing  rapid  variations  in  the  geomagnetic 
field.  Because  of  its  solar  origin,  this  noise  component  is  highly  correlated  with  the  1 1 -year 
sunspot  cycle  and  is  particularly  troublesome  during  peaks  in  sunspot  activity. 

In  addition  to  geomagnetic  noise,  magnetometer  measurements  made  near  the  ocean 
surface  are  subject  to  the  second  important  noise  component,  ocean  wave  noise.  The  ocean 
surface  is  a  dynamic  interface  separating  air  and  electrically  conductive  seawater.  Interaction 
between  the  earth’s  magnetic  field  and  this  moving  interface  induces  electric  currents  in  the 
conductive  medium  which  are  responsible  for  measurable  surface-generated  magnetic  fields. 
Bergin  has  computed  the  magnetic  field  attributable  to  surface  waves  generated  by  winds  of 
various  velocities.*  According  to  this  calculation,  a  steady  30-knot  wind  generates  a  surface 
wave  field  which  produces  over  1  nT  rms  of  magnetic  noise  at  100  meters  above  the  surface 
in  a  frequency  band  from  dc  to  1  Hz.  Considerable  effort  has  been  directed  at  characterizing 
this  noise  source  in  an  attempt  to  mitigate  its  presence  and  increase  the  probability  of 
detecting  magnetic  fields  produced  by  sources  below  the  ocean  surface. 

This  report  describes  a  scheme  for  eliminating  part  or  all  of  this  surface  wave¬ 
generated  noise  component  by  means  of  adaptive  noise  cancellation,  with  reference  data 
provided  by  a  signal  proportional  to  the  ocean  surface  displacement  measured  at  a  single 
point.  An  experiment  performed  at  the  NOSC  Oceanographic  Research  Tower  is  reported. 
Results  agree  with  Weaver’s  model  for  ocean  surface^  wave  noise  generation  and  demonstrate 
that  significant  suppression  of  this  noise  is  possible.” 


THEORY 

The  adaptive  noise  canceller  (ANC)  is  diagrammed  schematically  in  figure  1 .  A  signal 
s  combined  with  uncorrelated  noise  ng  forms  the  input  to  the  primary  channel  of  the  ANC. 
The  reference  channel  receives  a  second  noise  component  nj,  which  is  correlated  with  ng 
through  some  linear  physical  process  which  may  or  may  not  be  known.  Let  this  physical 
process  be  denoted  by  the  transfer  function  T(w).  Optimum  cancellation  occurs  when  the 


*nT  *  nanotesla  =  I0_<)weber/m-  *  10  ^  gauss  =  gamma  (7). 

*Naval  Research  Laboratory  Memorandum  Report  2843,  Magnetic  Variations  Caused  by  Wind  Waves, 
by  J.  M.  Bergin,  July  1974. 
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Weaver,  J.  T.,  Magnetic  Variations  Associated  with  Ocean  Waves  and  Swell,  Journal  of  Geophysical  Research, 
v  70,  p  1921-1929,  IS  Aprill965. 


3 


l 


1 


S 


.  .U-  ,  L, 


iI'Iim  ianVTSjij.-.  n. 


aril  —jAr au.-, 


adaptive  filter  takes  on  the  characteristics  of  a  Wiener  filter  with  a  transfer  function  given  by 


W« 


(1) 


where  SX(j  is  the  cross-power  spectrum  of  the  primary  and  reference  inputs  and  Sxx  is  the 
power  spectrum  of  the  reference.  If  the  reference  source  contains  no  signal  component  and 
is  uncorrupted  by  noise  which  is  uncorrelated  with  nQ,  then  it  can  be  shown  that  W(<a>) 
becomes  equal  to  T(oj).^  Adjustment  of  the  adaptive  filter  is  effected  by  a  gradient-seeking 
algorithm  which  attempts  to  minimize  the  total  power  of  the  error  signal 

e  =  s+n0-y  =  s+n0-h(t)*n  j  (t)  (2) 


where  h(t)  is  the  impulse  response  of  the  adaptive  filter.  For  our  purposes,  the  primary 
input  consists  of  a  combination  of  magnetic  noise  from  geomagnetic  and  wave  sources.  In 
a  real  application,  this  input  may  also  contain  a  signal  from  source?  of  interest  under  the 
ocean  surface.  Into  the  reference  input  passes  a  simultaneous  measurement  of  the  ocean 
surface  displacement  at  a  single  point  in  the  vicinity  of  the  magnetometer  sensor.  This 
iniormation  may  be  obtained  from  wavestaff  or  subsurface  pressure  measurements,  optical 
or  radar  range  measurements,  or  whatever  instrumentation  is  appropriate  to  the  platform  at 
hand.  Notice  that  magnetic  signals  arising  from  sources  under  the  surface  pass  through  the 
ANC  without  distortion  because  they  are  uncorrelated  with  the  surface  displacement 
reference. 

Although  the  adaptive  filter  finds  the  optimum  transfer  function  which  minimizes  c 
without  regard  to  the  physical  relationships  leading  to  W(oj),  it  is  instructive  to  derive  this 
expression  from  basic  physical  principles.  We  begin  by  focussing  on  a  simple,  sinusoidal 
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Figure  1 .  Adaptive  noise  canceller  (ANC)  schematic  diagram. 


^Windrow,  B„  et  al,  Adaptive  Noise  Cancellation:  Principles  md  Applications,  Proceedings  of  the  IEEE, 
v  63,  p  1692-1716,  December  197S. 
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ocean  surface.  Weaver  has  calculated  the  magnetic  Held  above  and  below  the  surface  that  is 
attributable  to  a  monochromatic  ocean  surface  gravity  wave  in  deep  water/  A  brief 
condensation  of  this  calculation  follows. 

Figure  2  shows  the  coordinate  system  applicable  to  this  problem  and  tab!?  1  lists  the 
notation,  initially,  let  us  assume  an  infinitely  deep  ccun  and  therefore  ignore  bottom 
effects.  The  induced  electric  field  in  seawater  is  V  X  (F  +  H).  The  earth’s  field  F  is 
enormous  compared  to  the  wave-generated  field  it,  so  that  the  source  term  above  is  just 
^  X  ?  to  a  good  approximation.  Maxwell’s  equations  for  the  field  vectors  may  then  be 
written  in  the  form 


?x£  =  -ait/at 

V  X  it  =  4»ro(f  +  ^  X  f ).  (3) 


To  solve  these  equations,  the  surface  velocity  field  V  must  be  specified.  Assuming 
that  the  ocean  medium  is  incompressible  and  irrotational,  and  that  the  amplitude  a  of  the 
wave  is  much  smaller  than  the  wavelength,  then  a  monochromatic  surface  gravity  wave  moving 
in  a  direction  0  from  the  x-axis  can  be  shown  to  have  a  velocity  given  by 


\?  *  -acj|"i(cos0i  +  sinflj)  +  k  exp  i  jcot-m(xcosd  4  ysinO)  exp(-mz) 


(4) 


where  m  *  ur/g.  This  surface  velocity  field  is  inserted  in  equation  3  and  solutions  are 
found  which  are  harmonic  in  time  and  displacement  for  regions  above  and  below  the  surface. 
Boundary  conditions  are  then  applied  to  these  solutions  which  require  that  the  normal  and 
tangential  magnetic  field  components  be  continuous  across  the  air-sea  interface.  Ultimately, 
this  procedure  leads  to  a  harmonically  varying  magnetic  field  whose  component  in  the 
direction  of  the  earth’s  field  (as  measured  by  a  total  field  magnetometer)  is  given  by 


f 


**  ,v-  ij r'lrlfcTkMM  fa1 1  rfA  flly  a 


(5) 


*>/ 


-iwaogF 

co 


S2  +  C2)e 


-sw2/g  i(wt-m(xcos0  +  ysin$>] 


where  wc  have  introduced  the  following  notation: 

S*  sin  I:  I  ■  dip  angle  of  the  earth's  field 

C  *  cos  1  cos  6  :  6  *  eastward  inclination  of  wave  propagation  from  the 
magnetic  meridian 

s  *  height  of  magnetometer  above  still  water  level. 

b f  is  Weaver’s  expression  for  the  magnetic  disturbance  attributable  to  a  perfectly 
sinusoidal  ocean  surface.  Of  course,  the  real  ocean  surface  is  a  complicated  superposition  of 
travelling  waves  of  variable  frequency,  amplitude,  phase,  and  direction  of  propagation. 
Equation  5  provides  the  magnetic  disturbance  associated  with  each  member  of  this  super¬ 
position.  By  means  of  a  complete  specification  of  the  surface  directional  spectrum,  one  rnay 
proceed  to  determine  the  total  magnetic  field  attributable  to  an  arbitrary  surface. 

Following  Pierson’s  treatment  of  the  ocean  surface  as  a  stochastic  system,  we  write 
the  surface  displacement  in  the  form 


0© 

P(x,y,t)  =  J 

— OO 


K  _ 

J  exp  i  [ wt-m(xcos£  +  ysin0)-rtto,0) )  '  v/[A(w,0))  “dddco 

-IT 


(6) 


where  t(cj  ,6)  is  a  phase  randomly  chosen  within  the  range  0  to  2ir  and  ( A(co,0)]  ~  is  the 
directional  spectrum  of  the  sea  surface.^  From  equation  5,  which  gives  the  magnetic  field 
component  parallel  to  the  earth’s  field  attributable  to  an  elemental  wave  of  amplitude  a.  we 
can  construct  a  similar  elemental  field  arising  from  wavelets  of  the  type  used  to  form 
p(x,y,t) : 


b t  -•  --agF  +  c:k^2/*  c* l tot~m(xe°s$  +  ysin0)-r(w.$>]  ’/[  A(w,0)] -dfldw 
w 


(7) 


Summing  all  such  contributions,  we  obtain  the  total  parallel  magnetic  field  component 
generated  by  p(x,y,t): 

rr 

f  -iirogF^T  +  £,2)e-sco2/g  ei[wt-mvxeos$  +  ysin#)-Tt'j,0)) 

J  00 


OO 

B/  (x.y.t)  =■  J 

-OO 


'  SlM  <x>,  #)l“d(5dcj  . 


(8) 


As  mentioned  previously,  the  reference  input  to  the  ANC  is  proportional  tc  'tx.y.t). 
The  primary  input  contains  the  noise  component  represented  by  B /  tx.y.i)  in  add.  .on  to 


^Pierson,  W.  J.  Jr.,  Wind  Generated  Gravity  Waves,  Advances  in  Geophysics,  v2,p  93-178,  1955. 
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other,  uncovrclated  noise  components,  Cancellation  of  B/  (x,y,t)  is  optimized  when  the 
adaptive  filter  assumes  the  Wiener  transfer  function  fdven  by  equation  I .  In  this  case,  the 
power  spectrum  of  tne  reference  input  is  just  the  power  spectrum  of  the  surface  displacement 
measurement  p(x,y,t),  which  may  be  obtained  by  first  forming  the  autocorrelation 
function 

oo  * 

<P(x,y,t)p*(x,y,t  +  X)>  =  J  J  lA(a>,$)l2e"lw*dwd$ 

-oo  -n  (9) 

where  we  have  assumed  that  {A( <*>,$))  represents  a  specially  homogeneous  and  temporally 
stationary  wave  field.  According  to  Pierson,  et  al,  (A(u>,6)]  can  generally  be  expressed 
ns  a  product  of  a  frequency-dependent  factor  and  a  direction-dependent  factor. 


lA(u>,#)]2  =  f(u>)cos20  ; 


■oo  <  oj  <  to 

-x/2  ^  9  ^  *12 


0  ;  otherwise 


(10) 


where  the  wind  direction  coincides  with  6  =  0.  Integrating  the  autocorrelation  function  over 
9,  using  the  general  form  foi  i  A(cu,0)]  given  in  10,  and  taking  the  Fourier  transform  yields 
the  power  spectrum  of  the  reference  input: 


Spp  =  ir/2f2(oo)  .  (11) 

Similarly,  the  primary  and  reference  cross-power  spectrum,  Spm,  may  be  calculated 
from  the  cross-correlation  function  of  the  surface  displacement  time  series,  equation  6,  and 
the  wave  generated  magnetic  field  time  series,  equation  8: 

<p(x,y.t)B^*(x’,y',t'+X)>  = 

7  f  i^(S2+C2)e-s"2/8  e-£’*  [A(oo,0)]Viw\j$dw 

■'  J  co 

-oo  -jr 

2 

— ^  CO  **  A. 

where  k-  r  = —  [(x  -  x')cos0 -t  (y  -  y')sin0]  . 

s  (12) 


The  factor  e“‘  ^  r  in  1 2  accounts  for  the  possibility  that  the  displacement  measurement 
may  not  be  taken  at  the  same  point  in  the  xy  plane  as  the  rr  .gnetic  measurement.  Fourier 
transformation  of  the  cross-correlation  function  gives  the  ctoss-power  spectrum  of  the 
reference  and  primary  inputs: 


$U,  S.  Navy  Hydrographic  Office  Publication  603,  Practical  Methods  for  Observing  and  Forecasting  Ocean 
Waves  by  Means  of  Wave  Spectra  and  Statistics,  by  W.  J.  Pierson  Jr.,  G.  Neumann,  and  R.  W.  Jones,  1955. 


(13) 


Som  *  (  (S2  ♦  C2)e"•<*,2/*  e"**1  ’  r  coa2#f2(<»>)d$ 

Pm  j  u 

-*n 

where  the  general  form  for  (A(u>,d))  2  given  in  10  has  been  used.  The  optimal  transfer 
function  is  then  given  by 

s  9,2 

e"*"2/*  f  (S2  ♦  C2)e-*^  *  ? cot2ld*  • 

SPP  “  -*/2  (14) 

Any  slow  variation  in  the  parameters  which  specify  this  transfer  function  will  be  sensed  by 
the  ANC,  and  the  adaptive  algorithm  will  adjust  the  transfer  function  of  the  adaptive  filter 
to  accommodate  these  variations.  This  feature  is  particularly  important  in  operational  settings 
where  the  magnetometer  sensor  heights  may  change  or  the  directional  characteristics  of  the 
wave  field  may  change. 


EXPERIMENT 

To  test  the  suitability  of  adaptive  noise  cancellation  techniques  to  the  reduction  of 
wave-generated  magnetic  noise,  an  experiment  was  conducted  at  the  NOSC  Oceanographic 
Research  Tower  located  about  1  mile  off  the  Southern  California  coast  near  San  Diego. 
Ordinarily,  sensitive  magnetic  measurements  of  this  type  are  hampered  by  stray  fields 
generated  by  steel  in  the  tower  structure.  At  this  facility,  however,  the  magnetometer 
sensor  may  be  mounted  approximately  20  meters  from  the  tower  on  a  nonmagnetic  boom 
made  of  aluminum  and  glass  fiber.  This  separation  provides  sufficient  magnetic  isolation 
from  the  tower  structure,  and  tests  with  sensitive  accelerometers  have  demonstrated  that 
boom  motion  is  negligible  in  all  but  the  severest  se.i  states.** 

Two  parallel  time  records  were  obtained.  One  time  record  monitored  the  total 
magnetic  field  as  measured  by  an  AN/ASQ-81  optically  pumped,  metastable,  helium 
magnetometer  mounted  7  meter,  above  the  surface.  This  instrument’s  maximum  sensitivity 
is  such  that  a  0. 1-nT  full-scale  signal  produces  a  15-volt  peak-to-peal;  output.  Under 
“quiet”  geomagnetic  conditions,  the  noise  floor  of  this  instrument  is  less  than  0.01  nT.  The 
widest  of  the  selectable  passbands  was  used,  giving  3-dB  points  at  0.04  Hz  and  2  Hz.  The 
second  lime  record  Wos  a  simultaneous  measurement  of  the  ocean  surface  displacement 
sensed  by  a  pressure  transducer  located  2  meters  below  the  still  water  level  and  10  meters 
north  of  the  magnetometer  sensor  head. 


**Podney,  W.,  and  R.  Sager,  Measurement  of  Fluctuating  Magnetic  Gradients  Originating  from  Oceanic 
Internal  Waves,  Science,  v  205,  p  1331-1382,  28  September  1979. 
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RESULTS 

The  Weaver  transfer  function  (equation  14)  was  derived  under  the  assumption  that 
bottom  effect*  could  be  ignored.  Such  conditions  hold  in  the  open  ocean.  However,  in  a 
coastal  environment,  wise  re  the  experimental  data  were  obtained,  the  presence  of  the  bottom 
should  be  included  in  any  theory  which  hopes  to  model  the  physical  process  responsible  for 
wave-generated  magnetic  Gelds. 

Modifications  to  Weaver’s  theory  which  account  for  finite  ocean  depths  have  been 
nisde  by  Woods. '  His  calculation  is  a  straightforward  extension  of  the  theory  outlined  fic»e 
and  need  not  be  described  in  dentil  except  to  say  that  an  additional  interface  between  the 
conductive  seawater  and  nonconductive  ocean  floor  must  be  introduced.  This  modified 
transfer  function  is  compared  here  to  the  experimental  data. 

Approximately  1  hour  of  ocean  tint  a  was  used  to  obtain  the  experimental  values  for 
the  transfer  function.  The  surface  displacement  itcord  was  divided  into  eight  consecutive, 
nonredundant  segments  of  400  seconds  in  length.  For  each  400-second  segment,  a 
windowed  Fourier  transform  was  computed.  The  eight  modified  periodograms  thus  obtained 
were  then  averaged  to  give  an  estimate  of  the  power  spectrum  of  the  surface  displacement. 

An  estimate  of  the  cross-power  spectrum  of  the  surface  displacement  and  the  magnetometer 
measurements  was  obtained  in  a  similar  fashion  and  the  resulting  transfer  function  was 
computed.  The  results  appear  in  figure  3,  where  the  solid  points  represent  the  experimental 


7  Defense  Research  Establishment  Pacific,  Victoria,  B.  C.  (Canada),  PNL  Laboratory  Note  65-6,  Magnetic 
Variations  Associated  with  Ocean  Waves  and  Swell  in  a  Shallow  Sea,  by  R.  S.  Woods,  July  1965. 


transfer  function.  Each  point  is  an  average  of  the  transfer  function  over  a  frequency  range 
which  spans  0.025  Hz  on  either  side  of  the  point. 

Plotted  on  the  same  axes  is  the  theoretical  transfer  function  of  Weaver  modified 
according  to  Woods  to  account  for  shallow  water.  For  the  purposes  of  this  calculation,  the 
following  parameters  were  used: 

I  =60° 

a  *  A  X  10**  ^  *  emu 
g  ■  980  cm/s^ 
s  *  700  cm 
F  *  50000  nT  . 

In  addition,  the  wave  vector  of  each  travelling  wave  component  of  the  surface 
displacement  was  assumed  to  be  refracted  normal  to  the  shoreline.  Thus,  6  -  90°  to  c 
good  approximation.  The  agreement  between  this  theory  and  the  experimental  data  is  well 
within  the  error  of  the  measurements. 

With  the  connection  between  the  sur face  displacement  and  the  associated  magnetic 
field  well  established,  the  next  step  was  to  determine  the  extent  to  which  this  magnetic  noise 
could  be  modeled  by  the  adaptive  filter  from  reference  data  provided  by  a  signal  proportional 
to  the  surface  displacement  measured  at  a  angle  point.  The  ocean  data  were  processed  with 
a  Rockwell  Universal  Adaptive  Filter,  which  is  a  digital  implementation  of  the  ANC  in 
figure  1.  A  full  description  of  the  LMS  adaptive  filter,  the  heart  of  the  ANC,  can  be  found 
in  the  extensive  literature  on  this  subject  (eg,  McCool  and  Widrow) Very  briefly,  the 
adaptive  filter  consists  of  a  tapped  delay  line  connected  to  the  inputs  of  an  adaptive  linear 
combiner.  The  adaptive  filter  is  shown  in  block  diagram  form  in  figure  4.  The  filter  output 
at  time  j,  denoted  by  yj.  is  the  discrete  convolution  of  the  input  veclcr 

xj  =  (Xj,Xj_j,  •  •  • ,  xj_n+j )  (15) 

where  the  components  of  this  vector  are  delayed  samples  of  the  input  signal  Xj,  and  the 
weight  vector 

Wr(Wij,W2j,  •  ,wnj)  (16) 

This  convolution  yields  the  filter  output 
n 

yj=  X  WU  xj-i+! 

i  =  1 

Thus  the  weight  vector  is  a  sampled  representation  of  the  impulse  response  of  the  adaptive 
filter.  The  real-time  length  of  the  tapped  delay  line  determines  the  maximum  impulse 
response  length  of  the  filter.  Within  this  limit,  the  filter  may  assume  whatever  impulse 
response  minimizes  the  error  cj.  Adaptation  of  the  weight  vector  is  achieved  with  the  LMS 


®NUC  TP  530,  Principles  and  Applications  of  Adaptive  Filters:  A  Tutorial  Review,  by  J.  M.  McCool  and 
B.  Windrow,  N.  r  b  ’977. 
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algorithm  which  updates  the  weights  according  to  the  prescription  given  by  the  following 
expression: 

Wj+I  -  V  •  <l7; 

H,  the  adaptive  feedback  constant,  influences  the  rate  at  which  the  filter  attains  the  optimal 
impulse  response. 

Several  filter  lengths  were  used.  The  filters  contained  from  64  to  512  taps  which 
represents,  at  a  sample  rate  of  5.8  Ilz,  a  real-time  filter  length  of  1 1  to  88  seconds.  Within 
this  range,  the  extent  of  wave  noise  cancellation  peaked  broadly  when  the  filter  consisted  of 
about  200  taps,  which  corresponds  to  a  real-time  length  of  35  seconds.  With  the  weights 
initially  zeroed,  approximately  5  minutes  was  required  for  the  weight  vectc  converge  to 
a  stationary  solution.  Figure  5  is  a  plot  of  the  rms  spectra  of  the  magnetometer  output  before 
(solid  line)  and  after  (dashed  line)  cancellation.  These  plots  are  taken  from  data  processed 
with  a  256-tap  filter  and  n  =  2“®.  Throughout  the  frequency  range  from  0  to  0.3  Hz,  where 
most  of  the  surface  wave  spectral  power  is  located,  the  level  of  noise  suppression  averaged 
about  10-15  dB.  At  the  main  swell  frequency,  about  0.08  Hz,  the  wave-generated  noise 
was  suppressed  by  about  20  dB.  Residual  noise  can  be  attributed  to  sources  other  than 
surface  waves,  with  the  bulk  of  the  unsuppressed  noise  coming  from  geomagnetic  sources. 
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Figure  5. 


CONCLUSION 

The  results  presented  here  demonstrate  that  ocean  surface  wave-genetated  magnetic 
noise  nay  be  significantly  suppressed  by  an  adaptive  noise  cancellation  scheme  which  uses 
reference  data  furnished  by  a  signal  proportional  to  the  surface  displacement.  Magnetic 
signals  originating  from  subsurface  sources  are  not  expected  to  be  affected  by  this  noise- 
reduction  technique  since  the  reference  data  are  related  solely  to  ocean  surface  behavior. 

Although  the  surface  wave  data  used  in  this  test  were  characteristic  of  coastal  waters, 
it  is  clear  from  the  Weaver  theory  that  similar  measurements  made  in  the  open  ocean  ought 
to  give  the  same  level  of  noise  cancelling  performance.  In  the  deep  ocean,  an  area  of  several 
thousand  square  miles  may  be  dominated  by  the  same  wind  conditions.  This  property  ensures 
that  the  adaptive  filter  transfer  function  will  be  time-invariant  on  a  scale  which  is  long 
compared  to  the  filter  adaptation  time. 
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Amplitude  of  ccfin  vave  at  the  surface 
Frequency  spectrum  of  the  sea  surface 
Directional  spectrum  of  the  sea  surface 
cos  1  cosfl 

Earth’s  magnetic  field 
Gravitational  acceleration 
Dip  angle 

t*>2/g 

sin  I 

Height  of  magnetometer  sensor  above  still  water  level 

Power  spectral  density  of  ocean  surface  displacement  measured  at  one  fixed  point 
Cross-power  spectrum  of  ocean  surface  displacement  and  magnetometer  output 
Eastward  inclination  of  direction  of  wave  propagation  from  magnetic  meridian 
Eastward  inclination  of  direction  of  wave  propagation  from  wind  direction 
Electrical  conductivity 

Magnetic  field  parallel  to  the  earth’s  field  attributable  to  a  sinusoidal  ocean  surface 


B/ 


Magnetic  field  parallel  to  the  earth’s  field  attributable  to  an  arbitrary  surface 


p  Surface  displacement 

T(u>)  Transfer  function  relating  surface  displacement  to  magnetic  disturbance 


Table  1 .  Notation  used  in  this  report. 
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